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By 
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Gennaro L. Goglia^ and Alok K. Venna^ 


1. INTRODUCTION 


Instability of viscous flow along a concave surface is in the form of 
counterrotating streamwise vortices commonly known as Goertler vortices. 
These three-dimensional vortices have been found to cause transition, indi- 
rectly, by interacting with other types of disturbances usually present in 
the flow. First, theoretical analysis of the problem was done by Goertler 
(ref. 1). A detailed survey of theoretical and experimental work done in 
this field was presented by El-Hady and Verma (ref. 2). However, little 
work has been done on the higher eigenstates of Goertler vortices. D'Arcy 
was the first to realize the presence of a second mode although his work was 
not published (ref. 3). Other efforts in this area include work by Aihara 
(ref. 4), and that of Herbert (ref. 5). Herbert presented the neutral sta- 
bility analysis of higher modes for an incompressible fluid, and came to the 
conclusion that the "most dangerous" wavelength of different modes are near- 
ly the same. The same conclusion is arrived at in the present work for 
compressible viscous fluids. Aihara and Sonoda (ref. 6) discussed the 
effect of pressure gradient on secondary instability of Goertler vortices. 

Experimentally, the presence of secondary instability of Goertler vor- 
tices was observed by Wortmann (ref. 7), Bippes (ref. 8), and Aihara (ref. 
9). Wortmann observed that secondary instabilities lead to distortion of 
the interface between two vortices and their oscillation in the spanwise 
direction. 
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The purpose oJ the present work is to investigete the effect of secon- 
dary instabilities of Goertler vortices with respect to their growth and 
amplifications. Mother aim of this work is to look into the effects of 
suction on higher modes of Goertler vortices. 

2. PROBLEM FORMULATION AND NIMERICAL PROCEDURE 

A comprehensive linear theory of the stability of compressible boundary 
layers over a curved wall was presented by El-Hady and Verma (refs. 2, 10 
and 11). A detailed formulation of the problem, and derivation of governing 
equations based on orthogonal curvilinear coordinates, is presented in 
reference 2. The governing equations of the problem are Navier Stokes equa- 
tions, continuity, energy and equations of state. (See ref. 2, equations 
9 thru 13). After appropriate nondimensionaliaation and linearisation, the 
six eight order system of equations are written as eight first order equa- 

tions, as follows; 
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®mn is a variable coefficient matrix. Nonzero elements of this 
matrix are given in reference 2. 

Similar suction parameter Yj used to study the effect of suction, is 
defined as 


(p*V*) 

Y - Ji-R (5) 

(p *U*) 

CC 


1 


1 


Equations (1), with boundary conditions (2) and (3), are solved using a 
variable stepsize integrator written by Scott and Watts (ref. 13). 


3. RESULTS AND DISCUSSION 
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Neutral stability curyes of 2nd and 3rd mode were studied for Mach 
nimbers from 0 to 4. Contours of constant growth rates were plotted for 
m* * 0.0 and 3.0 for two higher modes to study the growth of these 
higher eigenstates and to compare the amplitude ratio with that of the 
fundamental mode. The effect of suction has primarily been discussed for 
ut * 3. 0 (see section 3.5). 




3.1 Neutral Stability 

Figure 1 shows the plot of neutral stability curves (a » 0) for the 
2nd mode of disturbances for different Mach numbers. Figure 2 shows the 
same curves for 3rd eigenstates and it i5 obvious, by comparing the twa 
figures, that the family of curves move upward on the G, B plane for higher 
modes suggesting stabilizing effect of higher Mach numbers. The value of 

the Goertler number corresponding to the minimum of the curves is known as ^ 

• ^ 

critical Goertler number. This number represents the local stability 

limit. Figure 3 is a plot of critical Goertler number as Mach number for 
different modes. As the roach number is increased the critical Goertler 
number is not affected and remains constant until a point where it starts 
increasing. This is true for 2nd and 3rd mode; however, for the fundamental 
mode the curve increases continuously with the Mach number. This suggests a 
delayed stabilizing influence of conq>ressibility on higher modes. 



3.2 GRO'JTH RATES OF HIGHER EIGENSTATES 


Figures 4 and 5 show the contours of constant growth rates of 2nd and 
3rd mode for = 0.0. Figures 6 and 7 show the same contours for m^ “ 3.0. 

As can be inferred from the neutral stability curves (figs. 1 and 2), 
the family of curves for constant amplification rates move up in the G, 3 
plane for higher eigenstates. This indicates that the fundamental mode may 
be the most predominant of all existing modes as far as local ’Stability is 
concerned. 

Lines of maximum growth are shown in figures 4 to 7 by dash lines, and 
it can be thought of as the path along which vortices are mcJt likely to 
grow since a minimum of the curves corresponds to the most unstable wave 
number. 


3.3 EIGENFUNCTIONS OF HIGHER MODES 

Figures 8a to 8d give a comparison of the shape of 2nd eigenstates of 
disturbance quantities u, v, w and 6, respectively, at different Mach 
numbers for a wave number of 3 “ 0*6 and a growth rate of a “ 2. The 
corresponding Goertler numbers for m« * 0,1, 2, 3, and 4 are G ■ 14.57, 

13.757, 12.492, 12.499 and 14.683. The value of u, v, w and 0 are 
normalized with the maximum of u component for the corresponding Mach 
number. 

Figures 9(a-d) give a comparison of the shape of 3rd eigenstates for 

mac = 0,1, 2, 3 and 4. The flow conditions are the same as in figures 8(a- 

d). The corresponding Goertler numbers in this case are G = 31.01, 28.803, 

25.195, 23.699, and 25.771. It can be seen from figures 8 and 9 that the 

location of u , v . and 9 move away from the wall as the Mach 
max rain max 

number increases. This shows that the effect of compressibility on higher 
eigenstates is similar in nature to that of the fundamental mode, and the 
disturbances persist in a larger region as the Mach number increases. 

3.4 AMPLITUDE RATIO OF HIGHER MODES 

The critical Goertler number is not a useful quantity if transition 

prediction is required. Growth of vortices must be taken into account. The 
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total growth of vortices can be represented in terms of amplitude ratio, 
which can be defined in terms of G and a, as follows: 

3_ . A 1° lie (6) 

a 3 G G 
o o 

where 3 q and G^ are the amplitude and Goertler number corresponding 
to the neutral stability. 

Integration of equation (6) can be done along various paths. Figure 10 
shows the plot of amplitude ratio vs. Goertler number G for the first 
three eigenstates for mac * 0.0, It can be seen from this figure that 2nd 
and 3rd mode not only originate after 1st mode but also the rate of growth 
of amplitude (slope of curves) decreases for higher modes. This indicates 
again that the fundamental mode is the dominant mode at mac " 0,0. 

Figure 11 shows the same curve for ma * 3.0, and all previous obser- 
vations regarding figure 10 apply to it. In the above analysis the integra- 
tion was done along the maximum growth path until G “ 45. Experimentally 
the wavelength of the vortices have been found to be conserved along the 
length of curved surface (ref. 12). Lines of constant wavelengths can be 
shown to have tangent 3/2 on the G.B. plane and are defined by the 
nondimens ional parameter A, as follows: 


A = 




(7) 


This can be written in terras of G and 3 as 


A 


_G 

p3/2 


( 8 ) 


Amplitude ratios are calculated for constant A for different eigenstates 
by integrating equation (6) until G * 45 for mac “ 3.0. The result is 
shown in figure 12. 


It Is interesting to observe that the most unstable wavelength 
corresponding to the peak of the curves is the same for 2nd and 3rd mode, 

X 

and is approximately equal to 40. Ihe same quantitative observation was 
made by Herbert (ref. 5) for incompressible fluids. Large reduction in 
amplitude ratio is observed for higher modes. It can also be seen that the 
cut-off wavelength increases for higher modes. 

3.5 EFFECT OF SUCTION ON HIGHER EIGENSTATES 

Figure 13 shows the neutral stability curves for different values of 
suction parameter v, for 2nd and 3rd mode, for nhe * 3.0. The critical 
Goer tier number increases with suction in both higher modes. This can be 
seen more clearly in figure 14. However, for the fundamental mode the 
critical Goertler number first decreases and then increases as suction is 
increased. 

The overall effect of suction on the fundamental mode and higher modes 
still remains stabilizing, as can be seen from figure 11, where amplitude 
ratio has been plotted against G for 3.0 and v * -1.2. 

Figures 15(a-d) and 16(a-d) show the shape of eigenfunctions for three 
different values of v for 2nd and 3rd mode, respectively. It is clear 
from these figures that the effect of suction is in bringing the 
disturbances closer to the wall where viscous dissipation is higher and, 
thereby, stabilizing the flow. 


4. CONCLUSIONS 


Fundamental eigenstate of Goertler instability seems to dominate the 
flow, and higher eigenstates play a smaller role in controlling the growth 
of vortices of the different eigenstates, which are assumed to exist 
independently. 

Obviously, the actual phenomenon of Goertler instability is nonlinear 
ovnd would require an understanding of the interaction between different 
modes. However, the present linear theory does give us an insight into the 
problem. 
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The "most dangerous" wavelength foi fundamental and higher modes is 
nearly the same. 

Suction seems to have a locally stabilizing influence on higher modes 
contrary to the initial destabilizing influence observed on the fundamental 
mode. 

The overall effect of suction on higher modes is stabilizing* 
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Figure 1. Neutral stability curves of 2nd mode for different Mach 
numbers . 
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Figure 3. Effect of compressibility on local stability of higher modes 
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Figure 5. Contours of constant growth rates 
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Figure 9. Shape of eigenfunctions of 3rd mode for = 0 - 4. 
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Figure 10. Effect of higher inodes on amplitude ratio calculated along 
the locus of maximum growth rates at = 0.0. 
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Figure 11. Effect of higher modes on amplitude ratio calculated along 
the locus of maximum growth rates at M =3.0. Effect of 
suction (y = -1.2) on these modes. “ 
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Figure 13. Effect of suction on neutral stability curves 
for 2nd and 3rd mode for = 3.0. 
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Figure 14. Effect of suction 
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Figure 15. Shape of eigenfunctions of 2nd mode at - 3.0 for y 
-0.5 and -1.2. 
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